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ABSTRACT 



We present sensitive phase-referenced VLBI results on the radio continuum 
emission from the z — 1.87 luminous submiUimeter galaxy (SMG) GOODS 850- 
3. The observations were carried out at 1.4 GHz using the High Sensitivity 
Array (HSA). Our sensitive tapered VLBI image of GOODS 850-3 at 0.47" x 
0.34" (3.9 X 2.9 kpc) resolution shows a marginally resolved continuum structure 
with a peak flux density of 148 ± 38 /^Jy beam~^, and a total flux density of 
168 ± 73 /xJy, consistent with previous VLA and MERLIN measurements. The 
deconvolved size of the source is 0:'27(±0"12) x < 0:'23, or 2.3(dzl.O) x < 1.9 kpc^, 
and the derived intrinsic brightness temperature is > (5±2) x 10^ K. The radio 
continuum position of this galaxy coincides with a bright and extended near- 
infrared source that nearly disappears in the deep HST optical image, indicating 
a dusty source of nearly 9 kpc in diameter. No continuum emission is detected at 
the full VLBI resolution (13.2x7.2 mas, 111x61 pc), with a 4cr point source upper 
limit of 26 /xJy beam^^, or an upper limit to the intrinsic brightness temperature 
of 4.7 X 10^ K. The extent of the observed continuum source at 1.4 GHz and the 
derived brightness temperature limits are consistent with the radio emission (and 
thus presumably the far-infrared emission) being powered by a major starburst 
in GOODS 850-3, with a star formation rate of ~2500 yr~^. Moreover, the 
absence of any continuum emission at the full resolution of the VLBI observations 
indicates the lack of a compact radio AGN source in this 2; = 1.87 SMG. 

Subject headings: galaxies: individual (GOODS 850-3) — galaxies: active — 
galaxies: high-redshift — radio continuum: galaxies — techniques: interferomet- 
ric 
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INTRODUCTION 



The strength of the far-infrared (FIR) extragalactic background imphes that a large 
fr action of galaxy formation and black hole accretion activity is hidden by dust (see a review 



m 



Lagache. Puget. fc Dold l2005l ) . A large fraction (~ 30%) of the background has been 



resolved in the submillimeter into discrete hie 



1997; 


Hughes et al. 


1998; 


Bareer et al. 


1998) 



(e.g. 



Small. Ivison. &: Blain 



bright (> 5 mJy at 850 /im) submillimeter 
galaxies (SMGs) were found to be mostly at h igh redshifts of ^ ~ ^.-^^^^ they dominate 
the total star formation in this redshift range (jChapman et al.ll2003l . l2005l ). They have total 



infrared (IR) luminosities of 10^^ to > 10^^ Lq, corresponding to star formation rates of 
lO^-lO^ Mq jt-\ 

An outstanding question about bright SMGs is whether the FIR emission is powered by 
active galact ic nuclei (AGNs) or by starbursts. Th is has been studied at various wavelength s 



using X-ray (lAlexander et al 



2003 



20051). optical flBarger et al.lll999l : IChaDman et al.ll2005r) 



2008 



near-infrared (NIR) (ISwinbank et al.ll2004l ). and mid- infrared (MIR) data (iPope et al. 
Murphy et al.ll2009l ). At even longer wavelengths, radio and mm interferometers add another 
powerful way to investigate the existence of AGNs through morphology. A small sample of 
two SMGs were marginally resolved by the Submillimeter Array at sub-arcsec r esolutions, 
suggesting extended starbursts to be the main power source ( Younger et al.|[2008l ). By com- 
bining t he Multi- Elemen t Radi o Linked Interferometer (MERLIN) and the Very Large Array 
(VLA), I Chapman et al.l (120041 ) obtained high angular resolution (~ 0'.'3) 1.4 GHz images of 
12 SMGs. The majority (8 out of 12; 67%) of their sources were resolved by MERLIN-I-VLA, 
again suggesting extended starbursts. However, the resolution of MERLIN+VLA is insuf- 
ficient to probe radio emission from compact radio AGNs. Whether the remaining 33% 
of compact radio sources in the MERLIN+VLA sample are powered by compact nuclear 
starbursts or by AGNs thus remains to be tested. 

The high resolution of Very Long Baseline Interferometry (VLBI) observations permits 
a detailed look at the physical structures in the most distant cosmic sources. Also, sensitive 
VLBI continuum observations can be used to determine the nature of the energy source(s) 
in these galaxies at radio fre quencies. To date, several high redshift QSOs have been imaged 
at mi l liarcs econd resolution (IFrey et al.lll997l . l2003t iBeelen et al.ll2004l ; iMomjian et al.ll2004j . 
20051 . I2OO8I ). These are the highest resolution studies of such distant QSOs by far. Here 
we are expanding such work on high-z SMGs. In this paper, we present sensitive VLBI 
observations of the SMG GOODS 850-3. 



The source GOODS 850-3 (aka SMM J123618.3+621551, GN6. IChapman et al.l 12004 



-4- 



Pope et al. 2006 ). which has an optical redsh ift of 1 .86 j§ ( Chapman et al. 2005 ). is a lu- 
minous SMG with S'sso/xm = 7.72 ± 1.02 mJy (IWang et al.ll2004l ). This object is one of the 
stron gest radio sources among SM Gs, and it is designated as a compact source at 0.3" resolu- 
tion Jchapman et al.ll2004j . Eoosf ). The flux density of GOODS 850-3 at 1.4 GHz, measured 
with MERLIN-FVLA observations, is S'i.4ghz = 151 ± 11 /xJy. 

Throughout this paper, we assume a flat cosmological model with = 0.3, = 
0.7, and Hq = 70 km s~^ Mpc~^. In this model, and at the distance of GOODS 850-3, 
1 milliarcsecond (mas) corresponds to 8.4 pc. 



2. OBSERVATIONS AND DATA REDUCTION 

The VLBI observations of GOODS 850-3 were carried out at 1.4 GHz on 2006 December 
9 and 2007 February 13 using the High Sensitivity Array (HSA), which includes the Very 
Long Baseline Array (VLBA), the phased Very Large Array (VLA), and the Green Bank 
Telescope (GBT) of the NRAo|f|. Eight adjacent 8 MHz baseband channel pairs were used 
in the observations, both with right- and left-hand circular polarizations, and sampled at 
two bits. The data were correlated at the VLBA correlator in Socorro, NM, with a 2 s 
correlator integration time. The total observing time was 16 hr. In these observations, the 
shortest baseline is between the phased VLA and the VLBA antenna in Pie Town, NM, 
which is 52 km. This short-spacing limit filters out all spatial structure larger than about 
0''42. Table 1 summarizes the parameters of these observations. 

The observations employed nodding-style phase referencing using the calibrator J1229+6335 
(5'i.4GHz = 0.38 Jy) with a cycle time of 280 sec, 200 sec on the target source and 80 sec 
on the calibrator. The angular separation between the target and the phase calibrator is 
1.55°. A number of test cycles were also included to monitor the coherence of the phase 
referencing. These tests involved switching between two calibrators, the phase calibrator 
J1229-I-6335 and the phase-check calibrator J1219+6600 (S'i.4ghz = 0.06 Jy), using a similar 
cycle time to that used for the target source. The angular separation between these two 
calibrators is 2.6°. 



" Pope et al. 20081 derived a redshift of 2.00 ± 0.03 from Spitzer MIR spectroscopy with PAH features. 
Our observations in § 4.2 show that GOODS 850-3 is extremely optically faint and therefore the spectro- 
scopic observations of IChapman et al.l (|2005l ) might have been made on another galaxy. Nevertheless, the 
inconsistency between the two redshift measurements does not affect our analysis. 

•^The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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The accuracy of the phase cahbrator position is important in phase-referencing obser- 



vations (IWalkerlll999l ). as this determines the accuracy of the absolute position of the target 



source and any associated components. Phase referen cing, as used her e, is known to preserve 



absolute astrometric positions to better than ±0''01 (jFomalontlll999[ ) . 



Data reduction and analysis were performed using the Astronomical Image Processing 
System (AIPS) of the NRAO. After applying a priori flagging, amplitude calibration was 
performed using measurements of the antenna gain and system temperature for each station. 
Ionospheric corrections were applied using the AIPS task "TECOR" . The phase calibrator 
J1229+6335 was self-calibrated in both phase and amplitude and imaged in an iterative 
cycle. 

Images of the phase-check calibrator, J1219+6600, were deconvolved using two different 
approaches: (a) by applying the phase and the amplitude self-calibration solutions of the 
phase reference source J1229+6335 (Figure la), and (b) by self calibrating J1219+6600 in 
both phase and amplitude (Figure 16). The peak surface brightness ratio of the final images 
from the two approaches gives a measure of the effect of residual phase errors after phase 
referencing (i.e., 'the coherence' due to phase referencing). At all times, the coherence was 
found to be better than 93%. 

The self- calibration solutions of the phase calibrator, J1229-f6335, were applied on the 
target source, GOODS 850-3, which was then deconvolved and imaged at various spatial 
resolutions by tapering the visibility data. 



3. RESULTS & ANALYSIS 

Imaging the target source at the full resolution of the VLBI array, which is 13.2 x 7.2 mas 
(111 X 61 pc, PA=— 16°), achieved an rms noise level of 6.5 /iJy beam~^, but it did not reveal 
any continuum component in the field of GOODS 850-3. This indicates the absence of any 
compact radio continuum emission with flux densities of > 4cr or 26 yuJy beam~^, which in 
turn implies an upper limit to the intrinsic brightness temperature (corresponding to a rest 
frequency of ~ 4 GHz) of 4.7 x 10^ K for any compact radio source in GOODS 850-3. Our 
coherence tests during these observations using two VLBI calibrators show that the lack of 
a strong point source in GOODS 850-3 at the full resolution of the array cannot be due 
to the phase referencing procedure. The VLBI flux limit reported above is almost an order 
of ma gnitude lower than th e flux measured by the VLA+MERLIN (151 ± 11 /iJy) at 1.4 



GHz ( IChapman et al.l |2004[ ) . This immediately implies that more than 90% of the radio 



continuum emission in GOODS 850-3 is extended and not confined to a central AGN. 
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In the following, we assess whether our HSA observations can recover the flux seen 
by the VLA and MERLIN. To do so, we applied two dimensional Gaussian tapers on the 
visibility data with various values. However, the only image with a reliable ~ Aa, detection 
was obtained by applying a Gaussian taper falling to 30% at 0.5 MA in both the u- and v- 
directions. This gave a beam size of 0''47 x 0''34 (3.9 x 2.9 kpc^, P.A.=18°). The resulting 
image is shown in Figure 2. This image is practically made from the shortest baseline in our 
data set (phased VLA — Pie Town), but only after calibrating its visibilities using all the 
antennas in the array as described in §2. The rms noise level in this naturally weighted image 
is 38 yuJy beam"^. At this resolution, a Gaussian model fitting reveals a marginally resolved 
continuum source with a peak flux density of 148 ±38 yuJy beam~^ and a total flux density of 
168 ± 73 jUJy, which agrees w ell with the flux density measured with the VLA±MERLIN at 1.4 
GHz (jChapman et al.l 120041 ). Deconvolving the synthesized beam from the Gaussian fitting 
model of the source results in 0''27 ± 0''12 for the size of the major axis, and an upper limit 
of 0^'23 for the size of the minor axis, i.e., 0"27(±0"12)x < 0!'23, or 2.3(±1.0)x < 1.9 kpc^. 
The derived intrinsic brightness temperature limit is Tb > (5 ± 2) x 10^ K. 



DISCUSSION 



4.1. Radio Properties 



We have detected 1.4 GHz emission from GOODS 850-3 (corresponding to a rest-frame 
frequency of ~4 GHz) using the HSA. At a relatively low spatial resolution (0"47 x 0''34; 
Figure 2) the limit on the intrinsic brightness temperature value of the detected continuum 
source is > (5±2) x 10^ K. Its flux density at this resolu tion is consistent with that measured 
with the VLA±MERLIN Jchapman et al.lbooi [2005[ l. This implies that the radio contin- 
uum emission at 1.4 GHz is confined to the exte nt of the structure seen in the VLA±MERLIN 
image at 0''3 resolution ( IChapman et al.ll2004l ). or to the extent of the deconvolved size scale 
reported in this paper, which is 0!'27(±0!'12) x < 0'.'23, or 2.3(±1.0)x < 1.9 kpc^ 

At the full resolution of our array (13.2 x 7.2 mas), the radio emission from GOODS 850- 
3 is resolved out and does not show any single dominant source of very high brightness 
tempe rature (< 4.7 x 10^ K). This is in contrast to the results obtained by iMomjian et al. 
(I2OO4J ) on a sample of three high-z radio-loud quasars imaged with the VLBA, namely J1053- 
0016 {z = 4.29), J1235-0003 {z = 4.69), and J0913±5919 {z = 5.11). In each of these 2 > 4 
quasars, a radio-loud AGN dominates the emission at 1.4 GHz on a few mas size scale, with 
intrinsic brightness temperatures in excess of 10^ K. 



Condon et al, 



(Il99l[ ) have derived an empirical upper limit to the brightness tempera- 
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ture for nuclear starbursts. For a frequency value of 4 GHz (our rest frequency), the resulting 
limit on the intrinsic brightness temperature is Tb < 4.8 x 10^ K, while typical radio-loud 
AGNs have brightness temperatures exceeding this value by at least two orders of magni- 
tude. These authors also present a possible physical model for this limit involving a mixed 
non-thermal and thermal radio emitting (and absorbing) plasma, constrained by the radio- 
to-FIR correlation for star-forming galaxies. The measured intrinsic brightness temperature 
limit for the radio source in GOODS 850-3 is consistent with the empirical upper limit value 
for nuclear starbursts. 

For GOODS 850-3, the derived intrinsic brightness temperatures from our VLBI ob- 
servations are typical of starburst galaxies. However, when compared to local starburst 
powered Ultra-Luminous IR galax ies (ULIRGs ; Sanders fc Mirabel 1996), such as Arp 220, 



Mrk 273, and IRAS 17208-0014 fISmith et al.l Il998l : ICarilli fc Tavloii l2000l : iMomiian et al. 



20031 . 12OO6I : iLonsdale et al.ll2006l ). the radio continuum emission from GOODS 850-3 is an 
order of magnitude greater in luminosity, but also an order of magnitude larger in extent. 

In starburst dominated galaxies, the radio continuum is the sum of supernovae (SNe), 
supernova remnants, and re sidual relativist i c elec trons in the interstellar medium, as shown 
in the model presented by ICondon et al. (119911). However, d etecting individual SNe in 
GOODS 850-3 at z = 1.87 is unlikely. iLonsdale et al.l (|2006[ ) reported the detection of 
49 luminous radio SNe in the prototype ULIRG Arp 220 with flux densities between 0.053 
and 1.228 mJy and typical upper limits on their linear extent of less then 1 pc. At the 
distance of GOODS 850-3, the flux densities of such luminous radio SNe would be between 
(5 — 115) X 10~^ /iJy. These values are two to three orders of magnitude lower that the rms 
noise levels achieved in our VLBI observations. 

Our radio continuum results with the HSA clearly rule out a compact radio-loud AGN 
in GOODS 850-3 and show physical characteristics consistent with an extreme starburst. 
However, we cannot completely rule out that the radio emission in this source is from a 
kpc scale radio-jet structure with a very faint compact AGN component that contributes to 
less than 10% of the total radio continuum emission, and hence falls below our detection 
threshold at the full resolution of the array. Therefore, we further extend the discussion 
regarding the nature of the power source(s) in this SMG, and whether it hosts a radio-quiet 
AGN, by looking into the optical and IR properties of GOODS 850-3 at various bands, as 
described in the following section. 
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4.2. Multi-wavelength Properties 

Evidence for an intensive starburst in GOODS 850-3 is also seen at various other wave- 



lengths. First, ICarrilli &: YunI ( 119991 ) presented a method of using the local radio-FIR 



correlation and the radio and submillimeter flux densities to estimate redshifts of SMGs. 



Barger. Cowie. fc Richardsl (|2000[ ) derived a simple redshift formula for Arp 220-like Spectral 
Energy Distribution (SEDs): z + 1 = 0.98(5'85o/>S'i.4)°'^^. With this formula, the measured 
redshift z = 1.865, and the radio and submillimeter fluxes of GOODS 850-3, we find its 
radio and submillimeter SED to be consistent with Arp 220 within 20%. This implies that 
this SMG excellently follows the same radio-FIR correlation as Arp 220, and suggests that 
its FIR emission is powered by a starburst. 



GOODS 850-3 is among the first SMGs with an X-ray detection ([Alexander et al.ll2003l ). 
It is detected at soft X-rays in the 2Ms Chandra image of the Chandra Deep Field-North 
(CDFN), but it is not detected at hard X-rays. It has an X-ray luminosity of ~ 0.5 x 10^^ 
erg s~^. These authors suggested that this source is powered by star formation, but they did 
not rule out a highly obscured AGN with low X-ray luminosity. 



Pope et all (120081) obtained a MIR spectrum on GOODS 850-3 (GN06 in their paper) 
and decomposed the spectrum into PAH and "continuum" components. The continuum com- 
ponent contributes ~ 50% to its MIR luminosity, and such a component can be powered by 



a dust y AGN. The MIR spectral decomposition method was further refined by lMurphy et al. 



(120091 ). who noted that such MIR continuum fraction is a strict upper limit on the AGN 
contribution at this wavelength band. However, AGN dust components are warmer than 
those powered by starbursts and contribute mainly to the MIR. We therefore expect the 
AGN fraction in the total IR luminosity of GOODS 850-3 to be negligible. 

None of the above observational evidence, including our HSA observations, absolutely 
rule out the existence of an obscured AGN in GOODS 850-3. However, a kpc-scale starburst 
is consistent with all observations from the X-ray to the radio. 

The last piece of evidence comes from the optical and NIR morphologies. For the optical, 
we combined the latest version of the four (F435W, F606W, F775W, a nd F850LP) HST ACS 



images from The Great Observatories Origins Deep Survey (GOODS: iGiavalisco et al.ll2004| ) 
to form a deep "white" optical image for GOODS 850-3. For the NIR, we used three sets of 
data. First, Lihwai Lin et al. (in preparation) obtained nearly 30 hr of GOODS-N imaging 
data at J-band with the Wide-Field Infrared Camera (WIRCam) on the 3.6 m Canada- 
France-Hawaii Telescope (CFHT). We downloaded these J-band data from the public archive 
and reduced them. We also used the WIRCam on CFHT to obtain a deep i^^-band image 
with nearly 50 hr of integration. Lastly, we obtained a deep i^g-band image of the GOODS-N 
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with the Multi-Object Infrared Camera and Spect rograph (MOIRCS) on t h e 8.2 m Subaru 
Telescope. This MOIRCS imaging is described in IWang. Barger. fc Cowid ( l2009a| ). In this 
imaging, GOODS 850-3 received approximately 4 hr of integration. All the above J- and Ks- 
band images are still being deepened by various groups. All the data reduction was carried 
out by our group using SIMPLE Imaging and Mosaicking Pipeline (simpleE, I Wang et al.l 



m 



2009br). More details about the K.s-b and observa t ions an d data reduction can be found 



Barger. Cowie. fc Wand (120081 ) and IWang et al.l (l2009bl ) . The reduction of the J-band 



data is identical to that of the Ks-hsnd. Here we adopt the MOIRCS -ft'5-band image 
for morphological analyses because of its high angular resolution (FWHM ~ 0'.'4), and the 
WIRCam i^^-band image for photometry because of its greater depth. 

We present the ultradeep ACS white, WIRCam J, and MOIRCS Ks images of GOODS 850- 
3 in Figure 3, and the F435W to 8.0 /im photometry in Table 2. The ACS fluxes were 
measured with d = 0'.'8 apertures at the location of the Ks source. This is the maximum 
possible aperture size for the ACS fluxes to be free of contamination from nearby objects. 
The J-band image does not have sufficient resolution to separate GOODS 850-3 from nearby 
objects, and therefore we d id not attempt to measu re its J-band flux. The Ks-hand flux is 
measured with SExtractor ( Bertin fc Arnouti 1996 ) "AUTO" aperture, which approximates 
its total flux. In addition, Iwang et~aL J2009bh used a CLEAN-like method to construct 
Spitzer IRAC source catalogs based on the WIRCam Ks image. We include their IRAC 
fluxes for GOODS 850-3 in Table 2. We present the observed optical to IRAC SED in 
Figure 4. 

GOODS 850-3 is extreme ly faint in the optical, undetected by even the combination of 
the four ACS images (also see lChapman et al.ll2004l ). It shows some hint of flux at F850LP 
(see Table 2) but nothing appears visually in the F850LP image and the white image. It 
clearly shows up at J-band, however it is blended with a nearby galaxy. In the i^^^-band 
and Spitzer images, GOODS 850-3 is the brightest galaxy in its neighborhood. With the 
excellent, 0'.'4, resolution of Subaru at i^'^-band, the stellar component of GOODS 850-3 
is resolved. A simple Gaussian fit to the source in this Ks-hand image gives a FWHM of 
1'.'03 X 0'.'63, corresponding to 8.7 x 5.3 kpc^ at its redshift. The image shows an elongated 
morphology, which may be indicative of an edge on rotating disk, and a pronounced central 
region that coincides with the radio position, but does not dominate the total Ks luminosity. 
These are consistent with a nuclear starburst hosted by a massive galaxy. 

The above optical and near-IR observations imply a large galaxy of roughly 9 kpc in 
size that is entirely hidden by dust. This shows that dense molecular clouds are widely 



''also see 



http : //www . asiaa . sinica . edu . tw/ $\sim$whwcing/ idl/SIMPLE 
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distributed in this galaxy, despite the fact that only its nucleus is actively forming stars. 



4.3. Star Formation Rate and Surface Density 



The results obtained from various bands in the electromagnetic spectrum (radio to the 
X-ray) suggest that the dominant power source in the SMG GOODS 850-3 is a startburst. 
Therefore, in the following, we derive an estimate of its star formation rate and star formation 
surface density. 



The 850 yum flux density of GOO DS 850-3 implies a total IR luminosity of 1.5 x lO^^L© 



(LiR = 1.9 X lO^^S'sso/xm-^o mJy ^; iBlain et al.l |2002| ). Assuming that this IR luminos- 
ity is entirely powered by a starburst, then the associated star formation rate (SFR) is 
~ 2500 MfT. yr-| for a Salpeter initial mass function (SFR(M0 yr^^) = 1.7 x 1O"^°Lir(L0); 
KennicuttI 119981 ). An estimate of the SFR can also be calculated using the conversion 
factor between radio luminos ity and SFR (SFR(Mq yr'^) = 5.9 x lO'^s^i 4ghz(W Hz"!); 
Yun. Reddv. fc Condon! [2OO1I ) . which is derived using the local radio-FIR correlation. The 
radio luminosity of GOODS 850-3 is 3.7 x 10^^ W Hz"!, and the resulting SFR is ~ 2200 
yr~^. The consistency in the SFR values derived from the radio and the IR luminosities sim- 
ply confirms that this SMG follows the local radio-FIR correlation. Moreover, this suggests 
that the star formation activity is well confined within the radio source detected in our HSA 
observations. 

The size of the source seen in our 1.4 GHz image (Figure 2) compares well with other 
SMGs, but is co nsiderably smaller t han optically selected starbursting galaxies at high red- 
shift 



,see, e.g. 



Bouche et al.l 120071 ). From the size and the SFR of this SMG, we derive 
a star formation rate surface density of Ssfr ~ 700 Mq yr~^ kpc^^. This surface density 
is extremely high even after including the ~ 40% uncertainty in its size. For instance, the 
derived Esfr in GOODS 850-3 is a few times larger t han that seen in a sample of SMGs 



imaged at millimeter wavelengths (ITacconi et a. 



high end of SMGs (see Figure 3 of iBouche et al 



20061) . Such a high value is at the very 



20071 ). and is comparable to that measured 



i n the two extremely luminous SMGs resolved by high resolution submillimeter observations 
(lYounger et al.ll2008l ). Furthermore, the Ssfr in GOODS 850-3 is comparable to that mea- 
sured in the host galaxy of the z = 6.42 QSO SDSS J114816.64-I-525150.3 JWalter et al.l 



2OO9I ). Such high SsFR values are consistent with recent theoretical descriptions of (dust 
opacity) E ddington limited star f ormation of a radiation pressure-supperted starbursts on 
kpc scales ([Thompson et al.ll2005l ). 



Using the population synthesis model of iBruzual fc CharlotI (120031 ). the Hyperz package 



( iBolzonella. Miralles. &: Pelld l2000l ) package, and the ACS, Kg-, and IRAC photometry in 



Table 2, we find a stellar mass of 1.4 x 10^^ Mq in GOODS 850-3. The best-fit SED is shown 
in Figure 4. With a star formation rate of ~ 2500 Mq yr~^, the stellar mass can be doubled 
in just 56 Myr, roughly 1.6% of its Hubble time. The model also requires an extinction of 
Ay ~ 2.8 to explain the nondetections in the optical bands. We expect the extinction in the 
nuclear starburst component to be even larger. Furthermore, we derive a value of 18 Gyr~^ 
for the star- formation rate per unit stellar mass (i.e., specific star-formation rate; SSFR) in 
GOODS 850-3. This is an order of magnitude higher than that se en in normal star for ming 
galaxies at 2; ~ 2, which have SSFR values of about 2 — 3 Gyr^^ ( iPannella et al.ll2009[ ). 



In summary, GOODS 850-3 has a kpc-scale nuclear starburst with a star formation rate 
of ~ 2500 Mq yr~^ No evidence of an AGN is seen at any wavelength. An upper limit on its 
MIR AGN contribution of 50% is estimated by lPope et al.l ( 120081 ). but the AGN contribution 
to its total IR luminosity should be much lower. Moreover, our HSA observations reveal that 
more than 90% of the radio continuum emission in this source is extended and not confined 
to a central compact AGN. GOODS 850-3 has a stellar component of 1.4 x 10^^ Mq with a 
spatial extent of roughly 9 kpc. This huge and massive stellar component is entirely hidden 
by dust, suggesting an extremely rich molecular gas reservoir fueling the nuclear starburst. 
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Fig. 1. — Continuum images of the phase-check cahbrator J1219-I-6600 at 1.4 GHz: a) 
obtained by applying the phase and the amphtude self-cahbration solutions of the phase 
reference source J1229+6335, b) obtained by self calibrating J1219-I-6600 itself, in both 
phase and amplitude. The restoring beam size in both images is 17.1 x 6.4 mas in position 
angle —12°. The contour levels are at —3, 3, 6, 12, 24 times the rms noise level in the 
phase-referenced image (left), which is 1.71 mJy beam^^. The gray-scale range is indicated 
at the right side of each image. The reference point (0, 0) in both images is q;(J2000.0) = 
12''19°'35?7941, (5(J2000.0) = +66°00'3l'.'845. 
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Fig. 2.— Naturally weighted 1.4 GHz continuum image of GOODS 850-3 at 0.47" x 0.34" 
resolution (P. A. =18°). The peak flux density is 148 //Jy beam~^, and the contour levels 
are at —2, 2, 2.5, 3, 3.5 times the rms noise level, which is 38 /xJy beam^^. The gray-scale 
range is indicated at the right side of the image. The reference point (0, 0) is «(J2000.0) = 
12^36"^18?3317, 5(J2000.0) = +62°15'50'.'608. A two dimensional Gaussian taper falling to 
30% at 0.5 MA was applied. 
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Fig. 3.— Optical and near-IR images of GOODS 850-3. From left are the HST ACS "white" 
image (F435W+F606W+F775W+F850LP), the CFHT J-band image, and the Subaru Kg- 
band image. All panels have 10" sizes and north is up. Circles indicate the 1.4 GHz position 
and have 1" diameters. Note that GOODS 850-3 entirely disappears in the optical but 
becomes very bright at longer wavelengths. 



- 18 - 



100.00 r 

10.00 r 



^ 1.00^ 



— 1 1 1 1 — 1— 



—1 1 1 — I — — 



0.10 r 



0.01 r V 




10 



Wavelength {/^m) 



Fig. 4.— The SED of GOODS 850-3. Symbols show the observed data in Table 2. Filled 
squares are detections. Open squares are low significance flux measur ements. Downward 



triang les are 1 cr upper limits. The curve is the best-fit SED from the iBruzual fc Chariot 



(120031 ) models, which has a continuous starburst of 0.36 Gyr, extinction of Ay 
stellar mass of 1.4 x 1O^^M0. 
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Table 1. PARAMETERS OF THE VLBI OBSERVATIONS OF GOODS 850-3 



Parameters Values 



Observing Dates 2006 Dec. 9 & 2007 Feb. 13 

Total observing time (hr) 16 

Phase calibrator J1229+6335 

Phase-referencing cycle time (sec) 280 

Frequency (GHz) 1.4 

Total bandwidth (MHz) 64 



Table 2. OPTICAL TO MiD-lNFRARED PHOTOMETRIC DATA OF GOODS 850-3 



Wave Band Flux (/xJy) 



ACS F435W 
ACS F606W 
ACS F775W 
ACS F850LP 
WIRCam Ks 
IRAC 3.6 //m 
IRAC 4.5 //m 
IRAC 5.8 ^lm 
IRAC 8.0 //m 



0.012 ±0.010 
0.008 ±0.010 
0.030 ±0.015 
0.153 ±0.051 
4.71 ±0.14 

13.4 ±0.2 
17.9 ±0.2 
27.3 ± 0.6 

18.5 ±0.6 



